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Abstract 

In a previous work, we proposed an extended Nambu- Jona-Lasinio (NJL) model including heavy 
quark flavors. In this work, we will calculate strong and radiative decays of vector mesons in this 
extended NJL model, including light p, lo, K*, (j) and heavy D* , D*, B*, B*. 
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I. INTRODUCTION 

The Nambu-Jona-Lasinio (NJL) model [l|, l2], in its original form as a pre-QCD theory, 
was constructed of nucleons that interact via an effective two-body contact interaction. 
Late. *e .ode. wa. ..e.te,..e.ed as a *eo.. of ^ua.. degrees of freedom g Q • The .o. 
important feature of NJL model is the chiral symmetry of Lagrangian plus a chiral symmetry 
breaking ground state. The model was generalized to SU{3)f case of light quark flavors in 



refs 



.0 



a. 



On the other side, for heavy quark flavors the chiral symmetry no longer holds. However, 
new important symmetries such as the spin symmetry were discovered in heavy (Qg)-mesons 
[10|, which is a consequence of the order l/mg of spin-spin interaction in the effective quark 
potential ll|. In ref. 12|], The NJL model was generalized to include heavy flavors. Both 
the chiral symmetry in light meson sector and the spin symmetry in heavy meson sector 
were reproduced with the vector-current interaction. The bosonization technique was used 
there to obtain an effective Lagrangian of meson degrees of freedom. 

However as already shown in ref. |5|, vector-current interaction itself is not enough to 
reproduce the experimental masses of light vector mesons such as p, K* etc. Other chiral 
symmetrical interactions such as the axial- vector-current one, are needed to get satisfactory 
results for light meson sector. But these additional interactions do not obey the spin sym- 
metry in heavy meson sector since they will generate the incorrect spin-spin interaction that 



is not l/rriQ suppressed. In the above work [12|, the authors just introduced two coupling 



constants Gi and G2 for the ligth meson sector and another different coupling G3 for the 
heavy meson sector. 

In our previous work [l3|, we proposed a solution to extend the NJL model to comprise 
the heavy quark flavors. The NJL interactions were expanded with respect to l/ruf of con- 
stituent quark mass rrif just like the expansion in the heavy quark effective theory (HQET). 
Naturally the vector-current interaction is dominant while other interactions such as the 
typical axial- vector-current one should be l/rrif suppressed. We had performed numerical 
calculations for both the light and heavy meson sectors. The mass spectra flt the experimen- 
tal data quite well. The decay constants of heavy mesons were smaller than experimental 
values roughly by a factor of 2. 

The strong and radiative decays provide us important information about hadron struc- 



;ure. Experimentally, the decay widths of light vector mesons have been well measured 
14j-|l9| and so far some of decay widths or ratios of the charmed and bottomed heavy vector 

mesons were reported 20l-l22|. 

Generally speaking, it is a rigid test for any model to fit the experimental values of 

decay width or ratio. The most popular model for strong decay is the ^Pn model 23|, |24 |. 



This model has been applied to a great number of decay processes 



25 



. The radiative 



decays, mainly Ml transition which takes place when one of the constituent quark changes 



its spin and radiates one photon, has been studied in potential quark models 29|, |30| or from 



flavor symmetry 3l|]. For decays of h eavy mesons, abundant works have been done in the 



frameworks of chiral quark model 



30 



32|, potential model 



perturbation model [36|, and QCD sum rules |37|, l38|. The decays were also studied in NJL 



33 



34| . bag model |35|, chiral 



model 



39 



40| and from lattice QCD |4l|-|43|. 



In this work, we will calculate strong and radiative decays of vector mesons in the extended 
NJL model with heavy flavors, including light mesons p, u, K*, and heavy ones D*, D*, 
B\ B* 



II. MODEL AND FORMALISM 

In ref. [l3|, the Nambu-Jona-Lasinio model was generalized to deal with heavy quarks as 
well as light ones. The Lagrangian reads 



C = '4>{i$ - mo)^ + A, 



where 



A = Gy(^A^7/.V')' + 



mqUlqi 



[(^a^7m^)' + {i'Ki.i^i^f 



(1) 



(2) 



describes the four-point quark-quark interaction compatible with QCD chiral symmetry. 
Gy, of dimension (mass)~^, and the dimensionless h were parameters fixed in the spectral 
calculation. The second term on the right side in Eq. (|2]) appears as higher order correction 
expanded with respect to the constituent quark mass iriq similiar to the HQET expansion. 
We can rewrite Eq. ([2]) in a Fierz invariant form. For the light sector, one has 



a=i^Gvm)qf + {qil,\)qf] 



-7:{Gv + 



h 



rnqniqi 



)m)l,q? + {q>^)l^.l^qf 



(3) 



where AVs are the Uf{3) generators, with X'l = y |/ (where / is the 3x3 unit matrix) and 
the rest are Gell-Mann matrices in flavour space. For the heavy sector, one has 

A^ = iGviiQqy + (g^75A}g)(g^75A}Q)] 

- ^{Gv + -^^mi,q){qi'Q) + {qi,l,q){qYl,q)], (4) 

where still we have 

TrA.A,- = 25,,. (5) 

One can see that actually we only consider the higher order l/mqiriQ suppressed interaction 
in vector and axial-vector channels and so the important chiral symmetry breaking vaccum 
(the ground state) is unchanged. 

Using Bethe-Salpeter equation (BSE), we obtained the meson masses via the correspond- 
ing T-matrix where the mesons appear as the poles of the T-matrix. The meson-quark 
coupling constants were also obtained by further expanding the T-matrix around the meson 
poles. 

In this work, we will use the effective meson Lagrangian to calculate strong and radiative 
decays of vector mesons. The effective meson-quark coupling constants will be directly 
taken from our previous work. In the cases of pseudo-scalar meson and vector meson, the 
corresponding effective quark couplings read 

Kq = - 97rqi>il5'ri^ ' TT ^^7/^75^^^ " d^Tz, (6) 

Lpg = - gpqi}-i^,Ti} ■ p^. (7) 

For the decay of a vector meson (V) into two pseudo-scalars (P), one has 

r(V^ ^ PP) = -^ [ d(l)^^^\M{V -^ PP)\\ (8) 

2mv J 

where / dcj)^'^^ = J (2^)32^ (2'k)^^2E (27r)"^(5"^(g— fci — ^2) is the standard two-body phase-space- 
measure. In the rest frame of the decaying meson, the decay amplitude of the vector meson 
can be write as 

Miy ^ PP) = e% = -e ■ T, (9) 

where e^ is the polarized vector of V meson. Then we have 

4 



p + q 




p-q 



FIG. 1: The Feynman diagram corresponds to the strong decay process. 



The strong decay process of a vector meson is shown in Feynmann diagram Fig. [H where 



^ = (^, 0), / = ^ = (^, K), and 
of the constituting quarks. Using the Feynman rules, one can write down the expression for 



q = '""^'"^ = [^, U), L = -^ ^-^ = [ ' 2 ^ ^ ^c}-, ana mi, m2, m^ denote the constituent masses 



the decay amphtude directly. One finds 

d'^p 



iT^ = -Tr 



X 



(27r) 



iWvl 



-K92 



^A^- 



-i{9i + 



9i 



< — (^ — mi 



nil + ^^3 
i 



hYlrX' 



(11) 



+ f - m^ '"" ' 1712 + 7713'^' '" f + (^-m2 

For the reaction of a vector meson decays into a pseudo-scalar and a photon (7), V^ — )• P7, 
the decay width can be expressed as 

1 



viy -^ P-i) 



2m 



V 



i(2)|Al| 



where the decay amplitude should take the form 



iM{V -^ P-f) = ee^" {¥)€*" {j)T^ 



flV 



(12) 



(13) 



The Feynman diagrams of radiative decay are shown in Figj2j We can write down the 
radiative decay amplitude 



rpflU —J'j. 



d'^p 



rWvl 



^X^- 



-iQy 



-rrii f + t-mi 



>^A9p^9p ^"" ' 



mi + mi2 



)ii,y 



+ 4 - m2 



+ Tr I -^igvi^X""^ z(^p + ^p— ^^)^75A^ 



(2vr) 



r/ii 



r/ii + m,2 ' 



X ; iQY- ' 



/ - 7712 f + (^-m2 



k2 





FIG. 2: The Feynnian diagram corresponding to the radiative decay processes. 

In the rest frame of decaying meson, we only need the space components of the tensor T*-' 
and it can be written as 



q^ij _ Jjl'j'- 



VPy 



Then we have 



T{V -^ P-f) 



ak, 
3m^ 



C \rjn |2 

-l-VP-yl ; 



(14) 
(15) 



where a ^ 1/137 is the electromagnetic fine structure constant. 

To calculate the loop integrals, we apply the three-momentum cut-off regularization 
scheme to the integrals. First, we define some useful quantities 



Ep{m) =A/p2" 



771^ 



<:^i,2 = + g°±-Ep(mi), 
^^3,4 = -q^ ±Ep{m2), 
^5,6 = -f±Ek{mz). 

The cjj's emerge as poles when the integral with respect to p" is performed. After we integrate 
out p°, the amplitudes can always be represented as spatial integrals 



T 



d?p 



(27r): 



2,4,1 



E 



N 



PO=i^i 



47r2 



1 2,4,6 



p^dp I dt \^ 



N\ 



PO=^i 



'0 j~i —Uj^M-^j)' 

where N represents the numerator of integrand. The 2-dimensional integral will be per- 
formed numerically by Monte Carlo integration method using the vegas routine from gsl 
library. 



TABLE I: Meson-quark coupling constants. 



9n 


gx 


9D 


gos 


5B 


gBs 




4.25 


4.32 


4.71 


5.03 


5.92 


6.69 




gn 


gx 


go 


gos 


9b 


5Bs 




1.56 


1.61 


2.04 


2.09 


2.84 


3.11 




9p/u) 


ff</. 


gx* 


5D* 


5DJ 


ffB* 


9B* 


1.29 


1.38 


1.31 


1.64 


1.83 


2.51 


2.89 



III. NUMERICAL RESULTS 



In the previous work [13|, we had calculated the pseudo-scalar and vector mesons, light 
and heavy, consistently in an extended NJL model with interaction given by eq. ([2]). The 
input parameters were the current masses of light quarks and the constituent masses of 
heavy quarks, the two coupling constants and the 3-dimensional cutoff. Numerically, the 
parameters were set to 

m°/^ = 2.79 MeV, m° = 72.0 MeV, 

rUc = 1.62 GeV, nib = 4.94 Gev, (16) 

A = 0.8 GeV, Gv = 2.41, h = 0.65. 

Using above parameters we obtained the constituent masses of light quarks 

mu = md = 392 MeV, m, = 542 MeV. (17) 

The obtained meson-quark coupling constants, which we need to calculate the strong and 
radiative decays, are given in Table HI We will use the experimental meson masses given by 



Particle Date Group J44| 



In Table [TTl we show the results for the strong and radiative decays of light vector 
mesons. As we can see, our results are in qualitative agreement with the empirical val- 
ues. Nevertheless, quantitatively our results are systematically smaller than the empirical 



values by a factor of 2 or 3. The discrepancy always occurs in the NJL ca^ 



culation as the 



model lacks the quark confinement mechanism. In the potential model 45|], generally the 
masses of light vector mesons p or K* lay above the constituent quark mass thresholds and 
still they are bound states due to the linear confinement potential. In our calculation, the 



TABLE II: Strong and radiative decay widths for light vector mesons. 



Decay modes 




This work 


Bernard 
[39] 


Empirical 

[44] 


p —7- vrvr 


MeV 


68.5 


52.0 


149.1 ± 0.8 


p^ -^ 7r±7 


keV 


21.9 


60.1 


68 ±7 


p^ -^ TT^J 


keV 


43.9 


- 


89 ±12 


UJ ^)- TT-J 


keV 


866 


762 


764 ± 51 


(j) -^ K+K~ 


MeV 


1.28 


- 


2.08 


^^KlK^s 


MeV 


0.86 


- 


1.46 


K*^ -^ {Ktt)^ 


MeV 


20.9 


57.3 


50.7 ±0.9 


K*^ -^ K^-i 


keV 


13.5 


92.0 


50 ±0.5 


K*o ^ kO^ 


keV 


31.3 


- 


117 ±10 



constituent masses of light quarks are intentionally tune larger so that the mesons are still 
bound states under the constituent quark mass thresholds, even without the confinment. In 
another NJL calculation 39|, the smaller constituent quark masses were used and the p and 
K* vector meson was found as the resonant poles. Then they suggested to account for the 
discrepancy by introducing a renormalization factor of roughly 2 into the light vector meson 
field after have taken the higher order meson loops into consideration. In comparison, the 



numerical results from ref. 39J are also listed in Table [TTl As we know, the amplitudes of 
triangle FejTiman Diagrams heavily depend on the quarks masses when the meson masses 
are close to the mass threshold. Our numerical study shows that to fit the experimental 
decay width of p demands that 2mu should be very close to rup and then the numerical 
result turns to be unstable. We guess that the confinement mechanism is important here 
for the light vector mesons as it is critical to their formation. 

Table IIIII shows the strong and radiative decay widths of heavy vector mesons. Table IIVI 
exhibits the branching ratios for charmed vector mesons. It can be seen that our results 
agree with the experimental values. As the empirical data are not complete, here we also 
list some of other model calculation and lattice calculation in the table for comparison. In 
Table llllt our decay width of D*^ is a little lager than the empirical one. Numerically this 
can be corrected by changing rric slightly, about 5 MeV larger. In Table IIV[ our resulted 



TABLE III: Strong and radiative decay widths for heavy vector mesons (all in unit keV). 



Decay Modes 


This work 


Kamal 


Goity 


Empirical 






[46] 


[30] 


[20,21,47] 


D*^ -^ D^ttO 


39.7 


25.9 


28.8 




D*^ -^ D^TT^ 


84.4 


58.8 


64.6 




D*^ -^ Z)±7 


0.7 


1.7 


1.4 




D*^ -^ ah 


124.4 


86.4 


94.9 


96±22 


D*0 ^ D^TT^ 


46.5 


42.4 


41.6 




D*o _^ j)0^ 


19.4 


21.8 


32.0 




D*^ -^ ah 


65.9 


64.2 


73.6 


< 2.1 MeV 


Dl ^ Dsi 


0.09 


0.21 


0.32 


< 1.9 MeV 


B*^ -^ B^-f 


0.25 


— 


0.74 




B*o -^ B^-f 


0.22 


— 


0.23 




Bt ^ 5,7 


0.10 


- 


0.14 






TABLE IV: Branching 


ratios for charmed vector 


mesons. 


Decay Modes 


This work 


Kamal 


Goity 


Empirical 






[46] 


m 


m 


B,*± _^ B>±^o 


31.8 


30.0 


30.3 


30.7 ±0.5 


D*± ^ Z)0^± 


67.7 


68.0 


68.1 


67.7 ±0.5 


D*^ -^ D^j 


0.5 


2.0 


1.5 


1.6 ±0.5 


jj*0 ^ ^0^0 


70.6 


66.0 


56.5 


61 ±2.9 


£)*0 ^ ^0^ 


29.4 


34.0 


43.5 


38.1 ±2.9 



branching ratios also are in agreement with the experimental data. Here the numerical 
results are less sensitive to constituent quark masses than that of the light meson sector. 
We may expect that the calculation of strong and radiative decays for heavy mesons are 
more reliable as it is well known that for heavy mesons the confinement is less important 
than the one gluon exchange coulomb potential. 



IV. SUMMARY 

We have used the extended NJL model with heavy flavors [13| to calculate strong and 
radiative decays of vector mesons. It should be noted that no extra assumption and free 
parameter was introduced into our present calculation. A reasonable agreement to the 
experimental data is obtained. The results of light vector mesons may indicate that a more 
complex quark structure should be considered for vector meson due to the confinement which 
is lacked in NJL model. 
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